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Abstract 

We consider extended versions of the Hubbard model which contain addi- 
tional interactions between nearest neighbours. In this letter we show that a 
large class of these models has a superconducting ground state in arbitrary 
dimensions. In some special cases we are able to find the complete phase di- 
agram. The superconducting phase exist even for moderate repulsive values 
of the Hubbard interaction U. 
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In this letter we consider extensions of the Hubbard model. We shall show that they 
are superconducting. The idea to explain superconductivity in the framework of strongly 
correlated electrons was proposed in and has subsequently been studied in numerous 
publications. Most of these investigations used approximate or numerical methods and so 
only a few exact results are known. One of these exact results is an algebraic approach 
based on so-called 77-pairs @ . It allows for the construction of states exhibiting off-diagonal- 
long-range-order (ODLRO). This important concept has been developed in ||. ODLRO is 
an adequate definition of superconductivity in arbitrary dimensions since it implies also the 
Meissner effect and flux quantisation 

Let us first consider the Hamiltonian of the Hubbard model 

n (u) = -t E E (sW + clc ja ) + u ek-t - hfai - \) ■ (i) 

Here Cj a are canonical Fermi operators which describe electrons on a <i-dimensional lattice, 
i.e. {cj a , q t } = Sji5 aT and Cj a \ 0) = where | ) denotes the Fock vacuum. 

In (0) j runs through all L sites of the d- dimensional lattice, (jl) denotes nearest- 
neighbour sites and U is the Hubbard coupling. n$ a = Cj a Cj a denotes the number operator 
for electrons with spin a on site j and we write rij = + rij^. 

This letter consists of two parts. In the first part we shall add general nearest-neighbour 
interactions to the Hubbard Hamiltonian (HD and analyze the three following questions: 

• When does this new Hamiltonian commute with the ^-operators 

v = E c 3i c 3i » v j = E c n c h ■ ( 2 ) 
3 3 



• When will (77^ I ) be an exact eigenstate of this Hamiltonian? 

• When will (?7 T J | ) be a ground state of this Hamiltonian? 

We shall answer all these questions for the general multiparametric Hamiltonian. For 

example, | ipN ) = [V ) I ) wn l De a ground state in a 9-parametric region (subject to some 
inequality) . 

What makes the state | ip^ ) special is the fact that it has been shown to have ODLRO 
(and thus is superconducting), 

/ / 1 t t 1 / \ N N 

where we have also taken the thermodynamic limit (N, L — > 00 with N/L fixed). 

The results of the first part show that superconductivity based on r/-pairs is a rather 
typical phenomenon. To see the relation between ^-pairs and the usual order parameter of 
BCS theory we express the //^-operator in terms of electronic operators c k(7 in momentum 
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space , c] CT = A= J2k e ^ kc \ai which yields = Sfc cj^ct^. This is just the BCS order 
parameter. 

In the second part of this letter we will consider two special cases of the general Hamilto- 
nian discussed in the first part. Both of these models are exactly solvable in one dimension 
and have ?7-pairs in the ground state (and thus are superconducting) even for moderate 
repulsive values of U. These results are valid in any dimension. 

In order to argue that superconductivity based on //-pairing is a generic rather than an 
exotic phenomenon, we consider the Hamiltonian TC(U) = Ho(U) + TCi, where TLq{U) is the 
Hubbard Hamiltonian ([![), and 

Hl = X Y { C )a C lv + C \a C i<r){ n 3~v + U 1-<t) 
<jl>,<r 

+Vj2(n j -l)(n l -l) + J z j:S*S? 

<ji> <ji> 

+ J -f E (4s + 4$) + Y J2 (sV.Vu^t + M&m) 
<ji> <ji> 

+ P E (( n iT - \)( n n - \)( n i - !) + ("IT " \)( n H ~ \)( n J ~ !)) 

+Q Y Kt - \)i n n ~ h( n n ~ h( n n ~ «)> 

<ji> 

+/i Y n j + h E( n iT - n tt) ( 4 ) 

3 3 

where the SU(2) spin operators S" 1 are given by Sj = riji), Sj = c^c^ and Sj = cLc^. 

The first term in ([|) is known as the bond-charge interaction, the second one is the nearest- 
neighbour Coulomb interaction. In addition we included a XX Z-type spin interaction with 
exchange constants J xy and J z between nearest neighbour sites. The relevance of the pair- 
hopping term Y for high-temperature superconductivity has recently been discussed in ||. 
Apart from a chemical potential /i and a magnetic field h we also added a three- and four- 
particle density interaction P and Q, respectively. 

The Hamiltonian (f|) with t — X and zero magnetic field h = is the most general 
Hamiltonian that one can write down that is hermitian, symmetric under spinflip (a — > — a) 
and conserves the total number N a = J2j=i n ja of electrons with spin a and the number 
of doubly occupied sites. The last requirement is quite natural if we consider exact re- 
pairing ground states, since these are eigenstates of the number operator for doubly occupied 
sites N 2 = J2j n ji n ji- To include the Hubbard model itself we allow t ^ X and have for 
convenience also included a chemical potential \i and allowed for nonzero magnetic field h. 

The conditions under which T~t(U) commutes with can be read of from the identity 

[h(v), rf] = 2(t - x) y (4Ai + c )Ai) 
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+ (2V-Y)J2(v}(ni-l) + vl(n J -l)) 

<jl> 

+2P fa}(nj T - \){n n - ~)+Vi(n n - \){n n - -)) 

<ji> 

+2/iEi (5) 

3 

Therefore, the condition for ?7-symmetry is that t = X, 2V — Y = 0, P = and \x = 0. 
Using (^) one finds that (rf) N \ 0) is an eigenstate of Tl(U) if t = X and 21/ = Y. The 
corresponding eigenvalue is 

UL ZL / P Q\ / PZ\ 

where Z is the number of nearest neighbours of a lattice site. 

Next, we want to determine under which conditions the ^-pairing state is actually the 
ground state of the theory. For this, we use the following lower bound for the ground state 
energy. Let {i/j a } be an orthonormal basis of the Hilbert space, then the ground state energy 
E satisfies 

E > mm I (ip a \H(U)\ip a ) - £ \(ipp\H(U)\ip a )\\ . (7) 

V J 

This lower bound has previously been used by Ovchinnikov to analyze ferromagnetic 
ground states continuing the work of Strack and Vollhardt |TD[ . By requiring that this lower 
bound equals in (||) one can prove Jll[ that (rf) N \ ) is a ground state of H{U) if V < 
and 

^> m «(|P| + 2M + 2 +2 | t | + 2V ; + + (8) 

This shows that if t — X and 2V = Y < 0, one can for every value of the parameters in 7i 
always make the ground state superconducting, by making the Hubbard coupling sufficiently 
small. 

The results of [T2|] suggest that even if we relax the conditions t = X and 2V = Y, we 



still will have a superconducting ground state for sufficiently small U. This is much harder 
to prove, since the explicit form of the ground state is not known, but one might get some 
indications from doing perturbation theory in t — X or 2V — Y. 

One might object that negative U does not represent the real physical situation, since U 
represents the Coulomb repulsion between electrons and that should be positive. However, 
we shall see in a moment that in two special cases the following happens: as soon as U 
becomes larger than the bound (H), the ground state becomes (rf) N ~ n \ ip n ), where \ip n ) is 
a 2rz-electron state without doubly occupied sites. As U increases, so does n until n = N. 
At this value of U the theory will no longer be superconducting. It is therefore quite likely 
that superconductivity will persist for a somewhat larger range of values of U than given by 
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After these general considerations let us now look at two models more closely. For these 
models the full phase diagram in arbitrary dimensions can be obtained. Furthermore, in one 
dimension we are able to calculate the complete spectrum. From now on we will set t = 1 
for convenience. 

The first model is obtained by setting all parameters in (|]) equal to zero except for the 
bond-charge interaction X, 



H(X, [/) = -££ U) c c l(T + clc ja ) (1 - X {n h . a + n^)) + U £(n iT - \){n n - \] 



(9) 



For general values of X this model has been discussed extensively by Hirsch |T3[ . He argued 
that for large densities of electrons (low doping) the bond-charge interaction leads to an 
attractive effective interaction between the holes which may even create Cooper-pairs of 
holes. Indeed, using a BCS-type mean-field theory he found a superconducting phase for 
small hole concentrations. 

Unfortunately, up to now even in one dimension there are almost no exact results available 
although it has been shown that a simplified version of Hirsch's Hamiltonian in one dimension 
indeed has a strong tendency towards superconductivity Jl7| . 



In the following we consider the Hamiltonian for the special case X = 1. As mentioned 
above, at this point the number of doubly occupied sites is conserved. Therefore the Hub- 
bard interaction acts as chemical potential for the doubly occupied sites. The dynamics 
become quite simple since the local Hamiltonian permutes bosons (i.e. empty sites and dou- 
bly occupied sites) with fermions (i.e. singly occupied sites) on neighbouring sites but not 
bosons with bosons or fermions with fermions. In the sector with no double-occupations the 
Hamiltonian reduces to the well-known U = oo Hubbard model. 

In one dimension this Hamiltonian can be solved exactly [O] by generalising the method 



applied in [|lq , |T6|l to the U = oo Hubbard model. The Hilbert space is divided into certain 
subspaces in which the Hamiltonian can be mapped onto spinless fermions with twisted 
boundary conditions. The twisting angle depends on the subspace considered. In this way 
the complete spectrum can be obtained. The ground state energy for a system of Aq (single) 
electrons and N 2 doubly occupied sites at U = is given by 

sm(N-\7T /L) , . 

E (N 1 ,N 2 ) = -2 } 1 / } . 10 
sin (tt/L) 

Note that this energy is independent of N 2 . In general, the ground state is highly degenerate 
as is known for the U = oo Hubbard model. After including the Hubbard interaction we 
can minimize this energy for a given total number N = N± + 2N 2 of electrons. We find the 

phase diagram shown in fig. [1]. For U < —4 (sector I) the ^-pairing state \rf) \ 0) is a 
ground state in agreement with the inequality @. In sector II where —4 < U < U C (N) = 

— 4cos(iW/L) the ground state is of the form (rf} \U — oa,Ni) where \U = oo, Aq) 
is the ground state of the U = oo Hubbard model for Aq electrons. Both these states 
exhibit ODLRO showing the existence of a superconducting ground state even for moderately 
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positive values of U. For U > U C (N) the ground state is that of the U = oo Hubbard model 
(for densities D = N/L > 1 it is the state obtained by particle- hole symmetry). 

The phase diagram in higher dimensions looks quite similar although we can not give 

an explicit expression for U C (N). From (§) we know that m | 0) is a ground state for 
U < —2Z. The rest of the phase diagram can be obtained by an argumentation analogous 
to that in M. The properties necessary for this argumentation to hold are 

• 77-symmetry, 

• conservation of the number of doubly occupied sites, 

• for U — the ground state energy is independent of N 2 . 

Another interesting special case of the general Hamiltonian (|j) is the supersymmetric 
Hubbard model f|[7]]. The Hamiltonian of this model is given by 

HsH = H° + U |> iT - l -){n n - l -) (11) 
where H° = - and 

H li = c/tCjK 1 - n ji - n n) + c] t q t (i - n 3i ~ n n) 
+\( n 3 ~ l )( n i ~ 1) + c \Ai c n c i\ + cn c 3\ c \Ai 

+ Kt " ~\) + (mi ~ \)(nii - \) ■ (12) 

This Hamiltonian corresponds to the choice X = t = 1, V = —1/2, J xy = J z = 2, Y = — 1 
and P = Q = fi = h = 0in(^j. Note that H® t also contains a Coulomb interaction term 
—Z(riji — h)(riji — |). TC° commutes with r] and r] t as defined in (|J). Actually, the model 
has a larger set of symmetries that form the superalgebra U(2\2) B. 

The ground state phase diagram has been obtained in 0. It looks very similar to fig. [1|. 
The main difference is the occurence of the ground state \t — J) of the supersymmetric t — J 
model jTj| that replaces \U = 00) in fig. |1]. In addition, the sector I exists for all U < 
(in every dimension) as can also be seen from the inequality (§). A superconducting phase 
exist for moderate positive U in all dimensions. The additional interactions of ( P"2"|) lift the 
degeneracies encountered in the model @. A more detailed discussion of the phase diagram 
of the supersymmetric Hubbard model has been given in 0. 

In conclusion, in this letter we have shown that 77-pairs provide a simple and clear mech- 
anism of superconductivity. 
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FIGURES 

FIG. 1. Phase diagram for the Hamiltonian (9) in one dimension. D = N/L is the particle 
density and U the Hubbard interaction. Qualitatively the phase diagram in higher dimensions has 
the same form. 
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